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ABSTRACT 

A thermochemical cure model for predicting 
temperature and degree of cure profiles in curing 
phenolic parts was developed, validated and refined 
over several years. The model supports optimization 
of cure cycles and allows input of properties based 
upon the types of material and the process by which 
these materials are used to make nozzle components. 
The model has been refined to use sophisticated 
computer graphics to demonstrate the changes in 
temperature and degree of cure during the curing 
process. 

The effort discussed in the paper will be the 
conversion from an outdated solid modeling input 
program and SIN DA analysis code to an integrated 
solid modeling and analysis package (I-DEAS solid 
model and TMG). Also discussed will be the 
incorporation of updated material properties obtained 
during full scale curing tests into the cure models and 
the results for all the RSRM nozzle rings. 


INTRODUCTION 

Historically the model has been used to optimize 
cure cycles and to evaluate the effects of discrepant 
cure cycles on a part’s temperature and cure 
response. This model was utilized to a great extent 
during the Reusable Solid Rocket Motor (RSRM) 
throat pocketing erosion investigation efforts. The 
intent of the cure modeling effort was to determine if 
the current cure cycle for the RSRM throat ring 
would provide adequate cure within current 
specifications. Finite element models of the nozzle 
ring cross-sectional geometry were generated w ith 
the l-DEAS software package and then analyzed 
using the TMG thermal generator analysis computer 
code. Validation of the model with thermocouple 
data from full scale parts was possible due to this 
effort. Material properties in the model were also 
updated as a result of lessons learned during the 
nozzle erosion investigation. 


BACKGROUND 

The model was originally developed using the 
System Improved Numerical Differencing Analyzer, 
SINDA heat transfer program as a basis. 

Considerable effort was required to develop the user 
programmed subroutines to model other property 
changes which occur during cure. The model 
includes sub-routines for heat transfer total heat, 
thermochemical reactions and degree of cure. 

The exothermic heat generation rate, 0 CXO is 
expressed by the following: 

0 e xo = m resin AH™ 6a/6t 

Where 5a/5t is the phenolic resin or cure reaction 
rate, m resin is the resin mass, and AH™ is the heat of 
reaction. During the cure of phenolics a chemical 
reaction occurs causing an exotherm within the part. 
This exotherm is seen in the thermocouple data when 
the internal temperatures in the part surpass the 
temperature of the boundary. The reaction rate for a 
chemical reaction with an Arrhenius temperature 
dependence is as follows: 

da/dt = ke' &RT (1-a) 

where, 

da/dt = reaction rate, Ibm/hr 
k = frequency factor, lbm/hr 
E = activation energy, Btu-lb-mole 
R = gas constant (1.987 Btu/lb-mole deg R) 
T = temperature, degree R 
The values for resin kinetics, i.e. activation energy, 
frequency factor and heat of reaction are determined 
from Differential Scanning Calorimetry (DSC) 
testing. A 35% resin mass fraction was assumed for 
the carbon phenolic prepreg. The Activation Energy 
defines the threshold temperature at which the part 
begins to exotherm, the Frequency Factor determines 
the rate of reaction and the heat of reaction 
determines the magnitude of the reaction. 


In order to accurately predict the exothermic and 
viscous effects which occur during cure, it was 
necessary to use a two series reaction primarily for 
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the exothermic phenomena in parallel with another 
single reaction to enhance the viscosity model. 



1) Resin -> a -> (3 Series 

2) Resin -> 8 Parallel 

The equation for total heat dictates the amount of 
exotherm which occurs in the part once the threshold 
temperature has been reached. 

The equation for total heat is as follows: 

Qtouu = Q a * d a /dt + Qp * d p /dt 

S1NDA model grids were generated using SINGEN 
and the material orientation; cross-ply and with-ply 
properties were determined by hand. The new TMG 
code uses the 1-DEAS solid modeling to define finite 
elements and the material orientation input 
determines the material direction, see figure 2. The 
run times for analysis have been reduced from 30-45 
minutes to less than 1 minute. Figure 5 shows the 
comparison of the S1NDA and TMG results for the 
RSRM throat ring. 


OBJECTIVES OF RECENT EFFORT 

All nozzle ring models required updates to 
incorporate updated material properties obtained 
during UUEC efforts. Also all nozzle rings needed to 
be checked to ensure the current configuration is 
being used. First and second wrap cure models for 
all RSRM nozzle rings using, l-DEAS and the 
Thermal Model Generator node for l-DEAS were 
created, see figure I . 


FIGURE 1. RSRM Nozzle Rings 


The new model uses I-DEAS Solid Modeling; a three 
dimensional model of the part with geometry input, 
material properties and three dimensional mesh 
generator. Three dimensional mesh thermal model 
generator; TMG is a thermal modeling and analysis 
tool using a geometry-based finite difference 
modeling approach. Considerable effort was required 
to develop the user programmed subroutines to model 
property changes which occur during cure. The 
TMG model includes subroutines for thermochemical 
reactions and degree of cure. The model post process 
capability provides output of temperature and degree 
of cure versus time and position and displays the 
thermal output in movie style, see figures 3 and 4. 

FIGURE 2. FE Model and Material Orientation 

Carbon Cloth 
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FIG l Rt. J. Post process 



FIGURE 4. Post Process 



MODEL DEVELOPMENT AND 
VERIFICATION (MX 4926) 

Initialh values of conductivity tor cured carbon 
phenolic were used in the model, these values were 
adjusted m retleet prepreg properties by matching 
thermocouple data from emulated 504 ring panels 
instrumented and cured at MSEC. Resin kinetics and 
thermal ^onduciiv its 'acre also adjusted to provide 
better correlation h\ matching the thermocouple data. 
Comparison ol actual and predicted thermocouple 


data for each panel was done to ensure reliability of 
model. During the Nozzle Throat Ring erosion 
investigation a number of tull scale throat rings were 
instrumented prior to cure. This provided a unique 
opportunity to validate the cure model against full 
scale data. Figures 6 through 8 show some examples 
of the excellent correlation achieved between 
predicted and actual thermocouple data for the throat 
ring, forward nose ring and aft exit cone first wrap. 

Figures. Comparison of SINDA and TMG 

1 it Wr»p Throat ring 



FIGURE 6. Comparison of Predicted and Actual 
temperature 
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FIGURE 7. Comparison of Predicted and Actual 
temperature 


Forward Noaa ring 1 at Wrap 



FIGURE 8. Comparison of Predicted and Actual 
Temperature 


AH Eitt Con* 1 at Wrap 
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The most common cure model uses are assisting in 
cure cycle selection for new materials and part 
configurations, determining state of cure should an 
anomaly occur, providing temperature deviation of 
part should an anomaly occur and finally to study 
effects of material properties on curing such as the 
resin content and resin advancement. It’s important 
to understand that the cure model is used only as an 
aide in the disposition process and has never been 
used as the only decision maker. During D.R. 
(Discrepancy Report) disposition tag end and X-ray 
data is used to substantiate cure analysis results. The 
model has been used to help evaluate several full 
scale UUEC rings, 60K, MNASA insulation and 
recent MNASA cure discrepancies. 
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RESULTS AND CONCLUSIONS 


The sensible approach to upgrade the cure model for 
RSRM nozzle rings was to recreate all of the models 
in IDEAS based on current drawings and material 
properties. The linking capabilities within IDEAS 
allows automatic updating of existing models based 
on design changes when the Master Model is 
updated. 
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